Microscopic three-body forces and kaon condensation in cold neutrino-trapped matter 
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We investigate the composition and the equation of state of the kaon condensed phase in neutrino- 
free and neutrino-trapped star matter within the framework of the Brueckner-Hartree-Fock approach 
with three-body forces. We find that neutrino trapping shifts the onset density of kaon condensation 
to a larger baryon density, and reduces considerably the kaon abundance. As a consequence, when 
kaons are allowed, the equation of state of neutrino-trapped star matter becomes stiffer than the 
one of neutrino free matter. The effects of different three-body forces are compared and discussed. 
Neutrino trapping turns out to weaken the role played by the symmetry energy in determining the 
composition of stellar matter, and thus reduces the difference between the results obtained by using 
different three-body forces. 
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I. INTRODUCTION 

Within a newly born neutron star or proto-neutron 
star (PNS) neutrinos are temporarily trapped in a time 
scale of several tens of seconds, which would undoubt- 
edly contribute to the overall chemical equilibrium, and 
may change the structure and composition of the star. 
Neutrino trapping is expected to have a strong influence 
on the equation of state (EOS) of /^-stable neutron star 
matter |l|, and, as a consequence, affects the theoretical 
predictions of the maximum mass and the corresponding 
radius. The EOS of the neutrino trapped stellar matter is 
particularly important for the evolution of a newly born 
neutron star and the mechanism of black hole formation. 
The effects of neutrino trapping in stellar matter and its 
implications for astrophysicalphenomena were explored 
by several authors [j|llESlEE0- In Ref.Q, Prakash 
et al. investigated systematically the structure and com- 
position of PNSs by using various theoretical models in- 
cluding the schematic potential model [j| , the relativis- 
tic field approach based on the Walecka model [To| . 
and the effective chiral model E2- It was argued that 
a black hole would be most likely formed in a neutrino 
diffusion timescale of « 10 s, if the maximum mass of the 
hot neutrino-trapped star is about 1.5Mm in presence of 
negatively-charged hadrons. In Ref.0, 0, the presence 
of neutrinos was shown to delay the onset of hyperons, 
thus changing significantly the composition of star mat- 
ter, and making the EOS stiffer. In Ref. Q, the evolution 
of PNSs with quarks was explored. 

In the interior of a neutron star (NS) or PNS, the 
baryon density could be as high as several times the value 
of the nuclear matter saturation density, and thus the 
matter is expected to become exotic. Among the pos- 
sible exotic phases in dense nuclear matter, the kaon 
condensation is a subject of great interest in nuclear 



physics, hadronic physics and neutron star physics [Til 
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The presence of a kaon condensation 
may have important consequences for determining struc- 
ture [3 El El El, cooling rates H3, EH E3, H3| and evo- 
lution of neutron stars |24| . The phase transition from 
normal matter to the kaon condensed matter is also ex- 
pected to affect the transport properties and the glitch 
behavior of pulsars [2{| . Other strange phases could also 
appear in competition or coexisting with kaons; in par- 
ticular, hyperons are expected to occur at 2 — 3 times the 
saturation density pj, before the onset of kaon conden- 
sation. 

In Ref. [26| . the effects of different three-body forces 
on the neutron star maximum mass were investigated 
within the Brueckner-Hartree-Fock (BHF) approach. In 
Ref. , the kaon condensation in neutrino-free NS mat- 
ter was studied in the same theoretical approach. It was 
found that the composition of NS matter is sensitively 
dependent on the nuclear symmetry energy and the pres- 
ence of kaons. In the present paper we shall extend the 
previous work to the study of neutrino trapped matter 
at zero temperature, focussing particularly on the inter- 
play between the roles played by TBFs, neutrino trapping 
and kaon condensation. Effects due to finite tempera- 
ture, which are important for the internal com pos ition 
of protoneutron stars, and their evolution IT7L 123 . l29 j . 
will be discussed within the Brueckner-Hartree-Fock ap- 
proach, in a subsequent paper [30j. . 

The present paper is organized as follows. In Sect. 2 we 
review briefly the theoretical models adopted in our cal- 
culations, i.e., the BHF approach, and the chiral model 
for kaon-nucleon interaction. In Sect. 3 we discuss the 
stellar matter composition both with and without neu- 
trino trapping, and the resulting EOS. Our numerical 
results are presented in Sect. 4. Finally, a summary is 
given in Sect. 5. 
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II. THEORETICAL MODELS 
A. Brueckner-Bethe-Goldstone theory 

In the present work, we employ the Brueckner ap- 
proach for asymmetric nuclear matter |3ll l32| to calcu- 
late the baryonic contribution to the EOS of the stel- 
lar matter. The starting point of the BHF approach is 
the interaction G matrix which satisfies the Brueckner- 
Bethe-Goldstone (BBG) equation 31, 32] 

G{p : x pi uj) = v NN +v NN > — -— G{p,x p \u), 

(1) 

where lo is the starting energy, and x p = p p j p is the 
proton fraction, being p p , and p the proton and the total 
baryon density, respectively. Q is the Pauli operator, 
which prevents the two intermediate nucleons from being 
scattered into the states below the Fermi sea, and e(k) 
is the single particle energy given by e{k) = e(fc; p) = 
H 2 k 2 /(2m) + U(k; p). The single particle potential U(k) 
is calculated from the real part of the on-shell G-matrix 
and we adopt for it the so called continuous choice |33| , 
i.e. for any momentum k below and above the Fermi 
surface 



U{k;p) = Re (kk'\G[p;e{k) + e(k')]\kk') a , (2) 

k' <kf 

where the subscript "a" indicates antisymmetrization of 
the matrix element. Due to the occurrence of U(k) in 
Eqs. and @, the latter constitute a coupled system 
of equations that has to be solved in a self-consistent way. 
In the BHF approximation the energy per nucleon is 

J = l^ + Y E (kk'\G[p;e(k) + e(k')]\kk') a . (3) 

771 ? k,k'<k F 

Adopting the continuous choice for the single-particle po- 
tential, the two hole-line (BHF) truncation of the energy 
shift turns out to be a good approximation for the nu- 
clear EOS, since the results in this scheme are quite close 
to those obtained by including also the three hole-line 
contribution j34[. The basic input of the BBG equa- 
tion is the realistic baryon-baryon interaction, which is 
determined by reproducing the nucleon-nucleon scatter- 
ing phase shifts. The realistic nucleon-nucleon (NN) in- 
teraction vnn adopted in the present calculation is the 
Argonne wis two-body force |35|. Since nonrelativistic 
calculations, based on purely two-body interactions, fail 
to reproduce the correct saturation point of symmetric 
nuclear matter, three-body forces (TBFs) among nucle- 
ons are usually introduced §|6|. They turn out to be 
needed to correct this deficiency. In this work two mod- 
els have been used, i.e., the phenomenological Urbana 
model [33 , and a microscopic TBF constructed from the 
meson-exchange current approach |38| . In both cases, the 
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FIG. 1: The EoS is shown (panel (a)) for symmetric (lower 
curves) and pure neutron matter (upper curves). The sym- 
metry energy as a function of the nucleon density is displayed 
in panel (b). See text for details. 



BBG approach is applicable only if the TBF is reduced 
to an effective two body force, which has to be added 
to the bare two body interaction. This is done by av- 
eraging over the position of the third particle, assuming 
that the probability of having two particles at a given 
distance is given by the two-body correlation function. 
The resulting potential is a density dependent two-body 
force [3^|. The inclusion of the microscopic TBF in the 
calculation for asymmetric nuclear matter has been per- 
formed in Ref. |39j . 

In Fig.l, panel a), we show the equation of state for 
symmetric nuclear matter (lower curves) and pure neu- 
tron matter (upper curves). The solid lines represent 
the values obtained using the Urbana phenomenological 
TBF, whereas the dashed lines are the values obtained 
by adopting the microscopic TBF. We notice that the 
results of the two adopted TBFs show a similar behavior 
up to density p w 0.4 /to~ 3 , but they differ a lot in the 
high density range. In particular, the microscopic TBF 
turns out to be more repulsive than the Urbana model 
at high densities, and the discrepancy between the two 
predictions becomes increasingly large as the density in- 
creases. 

Nucleonic TBFs play an essential role in determining 
the chemical composition of a neutron star. In fact, as 
shown below, the baryon chemical potentials and the 
symmetry energy strongly depend on the choice of TBFs, 
for densities typical of those encountered in the neutron 
star cores. The proton and neutron chemical potentials 
can be derived from the energy density through the stan- 
dard thermodynamic relation 



P P ,n = (de/dp Pjn ) Pi , 



(4) 



with all remaining densities pi fixed. In asymmetric nu- 
clear matter, the calculation involves the knowledge of 
the energy per particle and its partial derivatives with 
respect to the total baryon density and proton fraction, 
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dE/A dE/A 
E/A + p^ x„ 



dp 



dx„ 



(5) 



dE/A , dE/A , . 

H P (p,x p ) = E/A + p^— + (l-x p )-^— (6) 



dp 



Actually, the energy per particle E/A may be expanded 
quadratically in the proton concentration x p about its 
value for symmetric matter (x p = }-)'■ 

|(A Xp) = |(P, ~) + (1 - 2x p ) 2 E sym + - 2x p ) 4 ) 

(7) 

Microscopic investigations 0, 0, have shown that the 
expansion up to the quadratic term is a good approxima- 
tion and the higher-order contributions is negligible. As 
a consequence, the difference of the neutron and proton 
chemical potentials is determined by the symmetry en- 
ergy in an explicit way 



ft,, 



fi p = 4(1 - 2x p )E sym (p) 



(8) 



which implies that the symmetry energy plays a decisive 
role in predicting the chemically equilibrated composi- 
tions of PNSs and NSs. 

In Fig.l, panel b), we show the symmetry energy cal- 
culated in the BHF approach as a function of the nucleon 
density. Its value at the normal nuclear matter density 
has been well determined to be 30 ± 4 MeV. Actually, 
the calculations performed with either the Urbana phe- 
nomenological TBF (solid line) or the microscopic TBF 
(dashed line) reproduce correctly this value. However, 
the high density behavior suffers a large uncertainty, and 
this reflects in the theoretical predictions of the neutron 
star composition. Moreover, we see that the symme- 
try energy predicted by using the microscopic TBF rises 
much faster at high densities than the one obtained by 
adopting the phenomenological TBF. 

As far as leptons are concerning, the energies and 
chemical potentials are obtained by solving the free Fermi 
gas model at zero temperature |40|. Once the baryonic 
and leptonic chemical potentials are known, one can pro- 
ceed to calculate the composition of the /3-stable stellar 
matter, and then the total energy per particle (E/A) to t 
and pressure P through the usual thermodynamical re- 
lation 



P = P 



d{E/A) tot 
dp 



(9) 



B. Kaon condensation 



In the present paper, we use the effective chiral the- 
ory proposed by Kaplan and Nelson and exten- 
sively investigated afterwards In this model, the 
SU(3) x SU(3) chiral Lagrange density is expressed as 



+ F Tri?7 P 7 5 [Afj, , B] + DTrB^ lb {A^B} 
+cTrM(U + [/ f ) + aiTrB(ZMZ + £m£)B 
+a 2 TrBB{£_M£,+^M^) (10) 
+a 3 TrBBTr{£MZ + ^M^). (11) 

The first four terms conserve the chiral symmetry, and 
the other terms break the chiral symmetry. Parameters 
in the symmetric part are / = 93 MeV, the pion de- 
cay constant, D = 0.81, and F = 0.44. The breaking 
strength is determined by the parameters 0.1,0,2,0,3, and 
c, as well as by the quark mass matrix M.. We adopt 
a i m s = —67 MeV and a 2 m s = 134 MeV, as deter- 
mined from the baryon mass splittings |42j . The con- 
stant c and the bare kaon mass are related by the Gell- 
Mann-Oakes-Renner relation m 2 K = 2cm s /f 2 . The pa- 
rameter 03 m s had remained largely uncertain for many 
years, due to our poor knowledge of the strangeness 
content of the proton and the kaon-nucleon sigma term 
"Ekn [HI EI [El [13 Fortunately, the large ambiguity in 
this parameter has been settled recently by Dong, Lagae 
and Liu 01 w ^ n small error based on the lattice calcu- 
lations. In the present calculations we adopt values of 
a 3 m s equal to —310, —222, and — 134 MeV, as done in 
Refs.[la.l23|. in order to investigate the sensitivity of our 
results to the variation of the proton strangeness con- 
tent. Our adopted central value a3m s = —222 MeV is 
very close to the value extracted from the lattice gauge 
calculations in Ref. 44], which is —231 MeV with error 
of less than 4%. The value of 03ms is essential in de- 
termining the onset density for kaon condensation, since 
it provides the attractive component of the kaon-nucleon 
interaction. 

In principle, the chiral Lagrangian should give the EoS 
of baryons and mesons, but so far only the kaon-nucleon 
interaction has been extracted from it. Therefore, we use 
the chiral Lagrangian only to extract the kaon-nucleon 
part of the interactions, and take the EoS for nuclear 
matter from the Brueckner many-body theory. This 
means that the energy density is a sum of three con- 
tributions: the kaon-nucleon interaction energy density 
(derived from C x ), the nucleon- nucleon energy density, 
and the energy density for leptons (e,/z): 

e = tKN + £NN + Clep- (12) 

From eq.Q, we get the nucleon- nucleon contribution as 



£nn{p,x p ) = p^(p,x p ), 



(13) 



whereas the energy density for leptons is the one for a 
free Fermi gas at zero temperature, and can be found in 
textbooks |40l |. 

Following exactly the standard prescription in 
Ref. |16(, we can get the kaon-nucleon energy density 
of the kaon condensed matter from applying the Baym 
theorem Bal. i.e. 



Trd^Ud^ + TrB^D^ - m B )B 



f 2 - 2/iio 2 r2 ■ 2 ® 

£kn = yMif Sin + 2m K f sin - 



(14) 
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+ p(2a\Xp + 2a 2 + Aa^)m s sin 2 — (15) 

where p and x p denote the nucleon number density and 
the proton fraction, respectively, and is the amplitude 
of the condensation. We remind the reader that in this 
work we neglect the effects of K° condensation, which 
could play a significant role at the high densities typ ical 
of neutron star cores, as it has been shown in Ref. |l8llT^ |. 



III. EQUATION OF STATE OF COLD 
CATALYZED STELLAR MATTER 

For stars in which the strongly interacting particles 
are only baryons, the chemical composition is determined 
by the requirements of charge neutrality and equilibrium 
under the weak processes 



Bi^B 2 + l+V u B 2 + l^B 1 



(16) 



where B\ and B 2 are baryons, I denotes a lepton, either 
an electron or a muon, and vi(yi) a (anti-)neutrino. Un- 
der the condition of neutrino escape, these two require- 
ments imply that the relations 



^qtXi + qiXi = 



Mi = biflr, 



qm 



(17) 



(18) 



are satisfied. In the above expression, x.- L — pi/ Pb rep- 
resents the baryon fraction for the species i, and ps the 
baryon density. The neutron chemical potential is de- 
noted by p n , whereas pi refers to the chemical potential 
of the baryon species i, bi to its baryon number and qi 
to its electric charge. The same notation holds true for 
those quantities with subscript I, i.e. leptons. Under 
condition when the neutrinos are trapped in the system, 
the beta equilibrium condition i|18|) is altered to 



Mi = & iMn - 9i(Mi - M*i) 



(19) 



where p Vl is the chemical potential of the neutrino vi. 

Because of the trapping, the numbers of leptons per 
baryon of each flavor of neutrino I — e, p, 



Yr 



(20) 



are conserved on dynamical time scales. Gravitational 
collapse calculations of the white-dwarf core of massive 
stars indicate that at the onset of trapping, the elec- 
tron lepton number Yj, e — x e + x Ve ~ 0.4, the precise 
value depending on the efficiency of electron capture re- 
actions during the initial collapse stage. Also, because 
no muons are present when neutrinos become trapped, 
the constraint Yl — x^ + x v = can be imposed. We 
fix Yl, at these values in our calculations for neutrino 
trapped matter. 
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FIG. 2: The particle populations are shown as a function 
of the nucleon density for neutrino free (left hand panels) 
and neutrino trapped matter (right hand panels). In the up- 
per (lower) panels results are displayed for the case when the 
phenomenological Urbana (microscopic) TBF is used. 
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FIG. 3: The particle populations are shown as a function of 
the baryon density for neutron star matter with kaons, for 
three different values of a^m B . In the upper (lower) panels 
results are displayed for the case when the phenomenological 
Urbana (microscopic) TBF is used. 



In the neutron star matter with kaons the chemical 
equilibrium can be reached through the following reac- 
tions 



p + I + vi , 



P + K~ 



I 



(21) 



where I denotes leptons, i.e., I = e,/i. One can deter- 
mine the ground state by minimizing the total energy 
density with respect to the condensate amplitude 9 keep- 
ing all densities fixed. This minimization together with 
the chemical equilibrium and charge neutrality conditions 
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leads to the following three coupled equations 16, 2 

cos0 = -—^[m 2 K f 2 + \up (2a 1 x p + 2a 2 +4 : a 3 )m s (22) 
/A 1 2 



a,m s =-310 MeV 



1 



and 



/i = fl e — /i„ c = [iK = Pn — ftp 



(23) 



(24) 



= 4(1 - 2x p )S{u) sec 2 - - 2aim s tan 2 - , (25) 



/ fi sin 9 + upo(l + x p ) sin — — x p upo 



Pe 



3tt 2 



= 0, 



(26) 



(27) 



where w = pb/ Pa is the baryon number density in units 
of pq. The last two equations are from the chemical equi- 
librium and charge neutrality conditions, respectively. 
For the neutrino- free case, the above equations recover 
the ones given in Ref.j^}. The EOS and the composi- 
tion of the kaon condensed phase in the chemically equi- 
librated neutron star matter can be obtained by solving 
the coupled equations (|23[1 . I|25l) . and 1271) . The critical 
density for kaon condensation is determined as the point 
above which a real solution for the coupled equations can 
be found. 

The stable configurations of a neutron star can be ob- 
tained from the well known hydrostatic equilibrium equa- 
tions of Tolman, Oppenheimer and Volkoff QUEUES for 
the pressure P and the enclosed mass m 



dP(r) 
dr 



Gm{r)£{r) 



1 



P(r) 
£(r) 



1 



4irr P(r) 
m(r) 



1 - 



2Gm(r) 



dm(r) 
dr 



4irr 2 £{r), 



(28) 



(29) 



once the equation of state P(£ ) is specified, being £ the 
total energy density (G is the gravitational constant). 
For a chosen central value of the energy density, the nu- 
merical integration of Eas. (|28l I29|) provides the mass- 
radius relation. For the description of the NS's crust, 
we have joined the hadronic equations of state above de- 
scribed with the ones by Negele and Vautherin ^| in the 
medium-density regime (0.001 fm~ 3 < p < 0.08 fm~ 3 ), 
and the ones by Feynman-Metropolis- Teller [43 and 
Baym-Pethick-Sutherland [f3(| for the outer crust (p < 
0.001 fm- 3 ). 



IV. RESULTS AND DISCUSSION 

In Fig. |21 we present the composition of the 
/3-equilibrated neutrino-free (left-hand panels) and 
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FIG. 4: Same as Fig. [^l but for neutrino trapped matter. 



neutrino-trapped npep, matter (right-hand panels). The 
upper panels display results obtained using the Argonnc 
vis two-body potential supplemented by the phenomeno- 
logical Urbana TBF, whereas in the lower panels we show 
the populations obtained when the microscopic TBFs are 
used. In the neutrino-free case, we notice that a stiffcr 
symmetry energy produces a larger proton concentration 
in the region p > 0.6 fm~ 3 , where the proton fraction 
x p w E 3 ym - I n neutrino-trapped matter, the lepton con- 
centration becomes sizeably higher because the electron 
chemical potential must keep at a high value, in order to 
fulfill the /^-equilibrium condition. This leads to a larger 
proton fraction as compared to the neutrino-free case, 
because of the charge neutrality. Therefore, neutrino- 
trapped matter is more symmetric than the neutrino-free 
one. 

We also notice that the differences in proton popula- 
tions, due to the different TBFs adopted, are smoothed 
out when neutrinos are trapped. This is readily under- 
stood as follows. In neutrino-free matter, the electron 
chemical potential is determined mainly by the density 
dependence of nuclear symmetry energy When neutrinos 
are trapped, they contribute to the chemical equilibrium, 
and, as a consequence, the electron chemical potential is 
determined simultaneously by the neutrino trapping ef- 
fect and by the symmetry energy. Because of the lepton 
number conservation, the major effect of trapping is to 
keep the electron concentration high, and this reduces 
considerably the role played by the symmetry energy in 
determining the star composition. 

Let us now discuss the case when kaons are present 
in (3 stable and neutrally charged matter. In Fig. [21 we 
show the particle fractions for neutrino free star matter, 
using alternatively the phenomcnological Urbana TBF 
(upper panels) and the microscopic one (lower panels). 
We have chosen three different values of the parameter 
a^rris, respectively a^rris = —310, —222, and — 134 MeV. 
We observe that the kaon threshold depends sensitively 
on the attraction term of the kaon-nucleon interaction. 
In the extreme case of aj,m s — —310 MeV, the pre- 
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dieted kaon onset density is the lowest, of the order of 
0.4 fm -3 , almost independently on the adopted nucle- 
onic TBF. This is due to the fact that, up to densities of 
the order of 0.6 fm -3 , both kinds of TBF produce a very 
similar symmetry energy. In addition, the role played 
by the kaon-nucleon interaction becomes more predom- 
inant over that by the symmetry energy for a stronger 
kaon-nucleon interaction term. For larger densities, the 
microscopic TBF gives a higher symmetry energy, and 
this leads to a slightly different kaon threshold when a 
weak attractive term is adopted. We also notice a dra- 
matic decrease of the lepton population and increase of 
the proton concentration as a consequence of the kaon 
appearance. This is due to the charge neutrality condi- 
tion, which has to be always fulfilled. We notice that 
the threshold densities for kaon condensation may be de- 
layed by the onset of hyperons, as found by other authors 

In FigQJthe particle population is plotted as a func- 
tion of the baryon density for the neutrino trapped case. 
We notice that, because of trapped neutrinos, the onset 
of kaon condensation is shifted to larger densities than 
in neutrino free matter. Even in that case, both TBF's 
produce very similar results, except for the less attrac- 
tive case, i.e. asiris = —134 MeV, where the onset for 
kaon condensation takes place at a lower density for the 
microscopic TBF. This is due to the fact that, at high 
density, the symmetry energy is larger, and this allows 
for a kaon condensed phase at baryon density lower than 
in the case with the Urbana TBF. As compared to the 
neutrino free case, the kaon abundance in the condensed 
phase is much smaller than in neutrino trapped matter, 
due to the larger electron population. In contrast to the 
normal phase, in the condensed phase the number of neu- 
trinos rises up as the baryon density increases, and this 
is because the presence of kaons leads to a decrease of 
the electron population as a function of the baryon den- 
sity. We notice that the values of the threshold densities 
for K~ condensation both in neutrino-free and neutrino- 
trapped matter are in very good agreement with those 
found in rcf. lj. However, the threshold densities are de- 
pendent both on the value of the temperature, and the 
attraction term of the kaon-nucleon interaction 0, Q . 

Once the relative particle concentrations are known, 
we can calculate the equation of state. This is shown in 
FigG3 both for the neutrino free case (upper panels), and 
the neutrino trapped one (lower panels). In the case with- 
out kaons, the equation of state of neutrino trapped mat- 
ter is slightly softer than the neutrino free one, because 
the loss in energy due to the reduction of the neutron- 
proton asymmetry as a consequence of the increase of 
the proton abundance exceeds the gain from the presence 
of neutrinos [l| . We observe that the kaon condensation 
produces a general softening of the equation of state with 
respect to the purely nucleonic case. The degree of soft- 
ening depends on the value of the parameter a^m s , i.e., 
the larger the value of |a3TO s | the stronger the softening 
is. In the case with kaon condensate, neutrino trapping 
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FIG. 5: The pressure for beta-equilibrated neutrino free (up- 
per panels) and neutrino trapped matter (lower panels) is 
shown as a function of the baryon density, for phenomenolog- 
ical (left-hand side) and microscopic TBF's (right-hand side). 
The thin solid lines represent the Maxwell construction. 
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FIG. 6: The gravitational mass is shown as a function of 
the radius (left-hand side) and the normalized central density 
(right hand side), for the phenomenological (upper panels) 
and microscopic TBF's (lower panels) in neutrino free matter. 
Circles represent unstable configurations. 



produces a stiffer equation of state due to the higher onset 
density of kaons and smaller kaon abundance, as shown 
in Fig0] It is evident from the figure that the soften- 
ing of the equation of state due to kaon condensate is 
much less pronounced in neutrino trapped matter than 
in neutrino free matter. This may lead a newly-formed, 
hot protoneutron star to metastability, i.e. a delayed col- 
lapse while cooling down, as discussed in ref.[Ell7|. 

In some cases, the onset of a kaon condensed phase pro- 
duces a negative compressibility in the equation of state. 



7 



Following Migdal [53|. we have performed a Maxwell 
construction to maintain a positive compressibility. We 
should note that the Maxwell construction is valid for 
substances with only one independent component. Neu- 
tron star matter has two independent components (the 
baryon and electric charge). Therefore, the Maxwell 
construction cannot satisfy Gibbs' criteria that pressure, 
temperature and all chemical potentials be common to 
both phases in equilibrium. A novel treatment of the 
phase transition in (5 stable matter has been proposed by 
Glendenning|54j for the hadron-quark phase transition, 
and has been subsequently extended to the study of the 
kaon condensed phase, treated as a first order phase tran- 
sition in neutron star matter [25j . The Gibbs construc- 
tion has strong consequences on the mechanical stability 
of neutron stars, because the pressure is a monotonically 
increasing function of the density. On the contrary, there 
is a mechanical instability for the Maxwell case that is 
initiated by the central densities for which the pressure 
remains constant. Such an unstable region is absent when 
the phase transition is treated using Gibbs' conditions. 
It turns out that the value of the maximum mass of neu- 
tron stars with kaon condensation is only slightly affected 
by the Gibbs construction[25j, which is still affected by 
many theoretical uncertainties Therefore we limit 

ourselves to the Maxwell construction. In Fig. 5, the re- 
gion of constant pressure (thin solid line) is comprised 
between two values of the baryon density, denoted by p\ 
and pi. If the magnitude of |ci3m s | is large, it may hap- 
pen that p\ < pq. As a rule, we consider this solution as 
unrealistic, and we do not show the corresponding curve 
in FigO This is the case when azm s = —310 MeV and 
neutrinos are trapped, that is to say that the kaon con- 
densation requires a very stiff symmetry energy if the 
kaon-nucleon interaction is strongly attractive. 

Once the EOS has been determined, the TOV equa- 
tions can be solved. The resulting gravitational mass is 
plotted for neutrino free matter in Fig|Hl as a function of 
both radius and central density, normalized with respect 
to the saturation value. In particular, microscopic TBF's 
produce stiffer equations of state, and therefore give rise 
to larger values of the maximum mass. We observe a 
decrease of the maximum mass when kaons appear, the 
exact value being dependent on the chosen TBF's and 
the parameter a^,m s . In the microscopic case, the inclu- 
sion of kaons decreases the maximum mass configuration 
down to about 2 M®. On the other hand, phenomeno- 
logical TBF's produce softer equations of state and, con- 
sequently, smaller neutron stars, whose maximum mass 
is about 1.7 M . In Fig|S] circles represent unstable con- 
figurations, originating from the Maxwell construction. 
Those configurations turn out to be stable if the Gibbs 
construction [2^, is applied to the transition from 
purely nuclear matter to the kaon condensed phase. 

In the purely nucleonic case, neutrino trapping gen- 
erally produces a softer equation of state, because beta- 
stable matter turns out to be more symmetric in neutrons 
and protons. As a consequence, the maximum mass be- 
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FIG. 7: The gravitational mass is shown as a function of 
the radius (left-hand side) and the normalized central density 
(right hand side), for the phenomenological (lower curves) 
and microscopic TBF's (upper curves) in neutrino trapped 
matter. 



comes smaller. This is displayed in Fig0 However, the 
appearance of kaons changes this general picture and the 
equation of state becomes stiffer. This is due to the fact 
that the K~ onset depends on the lepton chemical poten- 
tial, i.e., p e —p, Vc , which stays at larger values in neutrino- 
trapped matter than in the neutrino- free case, thus delay- 
ing the appearance of K~ to higher baryon density. The 
resulting maximum mass increases, as shown in Fig0 
If phenomenological TBF's are used, the value of the 
maximum mass stabilizes around 1.8 Mq, whereas micro- 
scopic TBF's give rise to protoneutron stars with mass 
slightly larger than 2 Mq. As compared to the neutrino 
free case (Fig|HJ), the reduction of the maximum mass 
due to the presence of kaons is smaller in the neutrino 
trapped case. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have investigated the properties of 
kaon condensed neutrino-free and neutrino-trapped mat- 
ter at zero temperature including the composition, the 
equation of state, and the radius-mass relation, in the 
framework of the BHF approach with both microscopic 
and phenomenological three-body forces. In particular, 
we have discussed the interplay among the effects of the 
different TBFs, kaon condensate, and neutrino trapping. 
It is found that the effect of the microscopic TBF re- 
sults in larger values of symmetry energy at high densities 
(pb > 0.4 fm~ 3 ) than that of the phenomenological one. 
This gives rise to a higher proton population in NS mat- 
ter without neutrinos and kaons, thus allowing for fast 
cooling through the direct URCA process. The contribu- 
tion of trapped neutrinos makes the electron concentra- 
tion keep high in /3-equilibrated star matter and weakens 
the role played by the symmetry energy on the predicted 
composition of the star. As a result, the composition of 
the star matter becomes less sensitive to the two different 
three-body forces considered when neutrinos are trapped. 
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In the star matter without kaons, neutrino trapping leads 
to a softening of the equation of state. However when 
kaons are allowed, neutrino trapping makes the equation 
of state stiffer as compared to the neutrino free star mat- 
ter, since the presence of trapped neutrinos shifts the 
onset of kaons to higher densities and reduces the kaon 
abundance. 

In both cases with and without kaons, the microscopic 
TBF leads to larger values of the maximum star masses 
than the phcnomcnological one. If only nuclcons and 
lcptons are allowed, the effect of neutrino trapping is to 
reduce the maximum star mass; while when kaons are 
allowed, trapping instead produces larger stars. As ex- 
pected, the effect of kaon condensate is generally to soften 
the equation of state and reduce the predicted maximum 
star mass. Neutrino trapping weakens the role of kaons 
by delaying their onset densities to higher values and by 
reducing the kaon abundance. As a result, the reduction 
of the maximum mass due to kaon condensate is less in 
the neutrino trapped case than in the neutrino free case. 

This general scenario could change if we include hyper- 
ons in our calculations. Several studies have been per- 
formed on hyperonic matter within the BHF approach 
|56l l57j. showing that E~ appear already at a very low 



value of the baryon density, i.e., pb ~ (2 — 3)/?o- There- 
fore, hyperon onset might be a process in competition 
with kaon condensati on pi, a s found in relativistic mean- 
field approaches |5lL l52t l58j|. We will present a more 
detailed study in a forthcoming publication. 
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